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ABSTRACT: Lymphoid enhancer binding factor 1 (LEF-1) is a member of the T-cell specific factor (TCP) family, which plays a key 
role in the development of breast endothelial cells. Moreover, LEF-1 gene has been identified as a candidate gene for teat number trait. 
In the present study, we detected two novel mutations (NC_010450.3:g. 99514A>G, 119846C>T) by DNA sequencing and polymerase 
chain reaction-restriction fragment length polymorphism in exon 4 and intron 9 of LEF-1 in Guanzhong Black, Hanjiang Black, Bamei 
and Large White pigs. Furthermore, we analyzed the association between the genetic variations with teat number trait in these breeds. 
The 9951 4 A>G mutation showed an extremely significant statistical relevance between different genotypes and teat number trait in 
Guanzhong (p<0.001) and Large White (p = 0.002), and significant relevance in Hanjiang (p = 0.017); the 119846C>T mutation 
suggested significant association in Guanzhong Black pigs (p = 0.042) and Large White pigs (p = 0.003). The individuals with "AG" or 
"GG" genotype displayed more teat numbers than those with "AA"; the individuals with "TC" or "CC" genotype showed more teat 
numbers than those with "TT". Our findings suggested that the 99514A>G and 119846C>T mutations of LEF-1 affected porcine teat 
number trait and could be used in breeding strategies to accelerate porcine teat number trait improvement of indigenous pigs breeds 
through molecular marker assisted selection. (Key Words: Pig, LEF-1 Gene, Expression Profile, Polymerase Chain Reaction- 
Restriction Fragment Length Polymorphism (PCR-RLFP), Teat Number, Haplotype) 



INTRODUCTION 

Porcine teat number is an important indicator of 
reproductive performance. Pig individuals vary in teat 
numbers and nipple development, which might be caused 
by a major gene and several minor genes (Wiesner et al., 
1994). However, the genetic model of the major gene and 
minor genes involved are still unknown. Therefore, it is a 
noteworthy and interesting task to study unreported genetic 
variations within candidate genes and evaluate their 
associations with teat number trait. 

Lymphoid enhancer binding factor 1 (LEF-1) is a 
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member of the T-cell specific factor (TCF) family, a small 
family including four members in vertebrates. As a key 
regulator of wingless signaling (Wnt), LEF-1 mediates P- 
catenin dependent transcription factors that bind to the 
promoter of several genes bearing the conserved sequence 
CTTTGT (Xie et al., 2005) and plays a critical role in the 
canonical Wnt signalling pathway required for 
morphogenesis of early mammary gland during 
embryogenesis. The LEF-l"'' mice lacked two pairs of 
placodes, and other pairs that formed were small and 
degenerated during mammary gland development (van 
Genderen et al., 1994; Boras-Granic et al., 2006). In 
addition, it has been recently shown that LEF-1 promotes 
oncogenic signaling in breast and other cancers (Arce et al., 
2006; Gebeshuber et al., 2007; Ravindranath et al., 2008). 
Misregulation of LEF-1 gene results in overactive Wnt 
signaling, driving LEF/TCFs to transform cells. Depletion 
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of LEF-1 causes down-regulation of cyclin D and over- 
expression of matrix metallopeptidase 7 in breast cancer 
(Wang et al., 2006) and breast cancer cell proliferation 
(Bucan et al., 2012). Therefore, LEF-1 is of particular 
interest to many groups for its significant effect on 
mammary gland development. 

As to its expansive role in mammary gland development, 
LEF-1 gene was identified as a candidate gene for teat 
number trait (Jonas et al., 2008; Tetzlaff et al., 2009). We 
tested the levels of LEF-1 mRNA in different tissues and 
the mammary gland at different developmental stages and 
found extremely high levels in the teat that had a tendency 
to gradually decrease during growth which suggested its 
enormous regulatory potential in mammary gland 
development. 

In previous studies, a 1351T>C transition and 1666A>C 

transversion in the 3' end of LEF-1 mRNA affecting porcine 
teat number trait in commercial populations were identified 
(Tetzlaff et al., 2009). However, in other pig breeds, 
whether the two mutants above exist and are there novel 
genetic variations are still unknown. Here, we aimed to 
identify novel genetic variations in the coding sequence 
region (CDs) and partly intron of porcine LEF-1 gene and 
analyze their associations with teat number trait of 
Guanzhong Black, Hanjiang Black, Bamei and Large White 
pigs. It was expected that the results of this study could 
provide some information for enriching pig genetic 
resources and accelerate the selection and breeding through 
marker assisted selection (MAS) to improve teat number 
trait in Chinese native pig breeds. 

MATERIALS AND METHODS 

Ethics statement 

Experimental animals were meticulously fed and the 
experimental procedure was in agreement with the Animal 
Care Conamission of the College of Animal Science and 
Technology, Northwest A&F University, China. 

RNA extraction and quantitative reverse transcriptase - 
polymerase chain reaction 

To survey expression of porcine LEF-1 gene in different 
tissues, glands and organs, quantitative reverse transcriptase 
polymerase chain reaction (qRT-PCR) was performed. Total 
RNA was isolated from teat, heart, liver, subcutaneous fat, 
spleen, kidney, lung and longissimus muscle tissues of 3- 
day-old Guanzhong Black pigs, teat and teat tissues of 45- 
day-old Guanzhong Black pigs as well as heart, liver, spleen, 
lymph, thyroid and fat tissues of 180-day-old Guanzhong 
Black pigs using TRIzol method. The integrity of RNA was 
checked by electrophoresis on 1.0% agarose gels with 
ethidium bromide (EB) staining. Three micrograms of total 



RNA from each category as well as control was reverse 
transcribed to obtain cDNA using moloney murine 
leukemia virus reverse transcriptase (Invitrogen) and 
ohgo(dT)18 primer following the manufactvffer's protocol. 
The primers for LEF-1 were 5'-AAATAAAGTGCC 
CGTGGTGC-3' (forward) and 5'-GAGAAAAGTGCT 
CGTCGCTGTA-3' (reverse), and P-actin (reference gene) 
primers were 5'-CGTGAAAAGATGACCCAGATCA-3' 
(forward) and 5'-CACAGCCTGGATGGCTACGT-3' 
(reverse). Each individual sample was run in triplicate wells 
and the amplification effects were about 95%. PCR 
amplification cycles were performed using iQ Multicolor 
Real-Time PCR Detection System (Bio-RAD, Hercules, CA, 
USA) and SYBR Premix Ex Taq II kit (TaKaRa, Dahan, 
Liaoning, China). The reactions were initially denatured at 
95°C for 3 min followed by 50 cycles of 95°C for 10 s, 
60°C for 15 s, 72°C for 35 s. The melting curve analysis 
was performed after amplification to verify the accuracy of 
each amplicion. The density of SYBR green I and 
determine the threshold cycle (Ct) value were analyzed by 
iQ Optical System Software 2.1. The change of transcript 
abundance of all tested genes was calculated using 2"^'^*'' 
method (Livak et al., 2001). All data were expressed as 
mean+SE. Accordingly, LEF-1 mRNA levels of all the 
tested tissues were given by the formula F = 2 AACt = 

(Ct, LEF-1 ~ Ct, p-actin)some tissue ~ (Ct, LEF-1 ~ Ct, p-actin)heart ^nd 

reported as fold change in abundance relative to heart 
LEF-1 level. 

DNA samples 

Three hundred and sixty-one ear samples were obtained 
from four pig breeds including Guanzhong Black pigs (GZ, 
N = 198, Hanjiang Black pigs (HJ, N = 34), Bamei pigs 
(BM, N = 30) and Large White pigs (LW, N = 99). Genomic 
DNAs were extracted from ear tissues using salt-chloroform 
extraction protocol (Sambrook et al., 2001) and stored at 
-80°C. The records of teat number trait from all subjects 
were collected for statistical analysis. 

Primer design and polymerase chain reaction 
amplification 

Based on GenBank sequence No. NC_010450.3, 13 
pairs of PCR primers (Table 1) were designed to amplify 
exon 1-12 and partial intron 1-11 of the pig LEF-1 gene. 
The 50 |iL volume contained: 400 ng genomic DNA, 2.5 
|iM of each primer, 25 |iL MIX (including MgCl2, dNTPs, 
Taq DNA polymerase and Bromophenol blue). The cycling 
protocol was 8 min at 95°C, 32 cycles of 94°C for 30 s, 
annealing (°C) for 30 s as indicated in Table 1 and 72°C for 
30 s, with a final extension at 72°C for 10 min. After 
amplifying the pig genomic DNA pool (2 (xL from every 
individual) (Sham et al., 2002), PCR products were 
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Table 1. Primer information for PCR amplification of porcine LEF-1 gene 


Loci 


1^1 illlCl SCL[LlCllCCa 


Tm 


OlZjCh allU lOCtlLlOU 


GenBank 


Detection 




Acc. No. 


methods 


PI 


F:CCAGCGGAGCGGAGATTACA 


57.0 


303 bp, including 


NC_010450.3 


DNApool 




R:AGCCTTCCCTCCCTGCTCTT 




5'UTR,El,partiaHVSl 


(nt32-334) 


sequencing 


P2 


F:TTCTTTCTTTTGTGGGTGGG 


60.0 


136 bp, including 


NC_010450.3 


DNApool 




R:AACGATAACATTGGGAGTGGTA 




partial IVSl; E2; partial IVS2 


(nt2679-2715) 


sequencing 


P3 


F:TTGTTCTTTGCATTCTCTCAGG 


59.0 


303 bp, including 


NC_010450.3 


DNApool 




R: GTC ATTCCTC ATAGC AGGTCTTA 




partial IVS2; E3; partly IVS3 


(nt4139-4442) 


sequencing 


P4 


F:TGATTCTTGTTACTGCTGATGC 


60.0 


344 bp, including 


NC_010450.3 


DNApool 




R:ACCACCCACCCAATGACA 




partly IVS3; E4; partial IVS4 


(nt81993-82337) 


sequencing; 












PCR-RFLP 


P5 


F:GAATGTCCTGACTCTGCCTCT 


59.0 


256 bp, including 


NC_010450.3 


DNApool 




R:TGCTTGTTGATGTGAGGTTTC 




partial IVS4; E5; partial IVSS 


(nt86310-86536) 


sequencing 


P6 


F:GAAGGGAAGCAGGCAATG 


59.0 


286 bp, including 


NC_010450.3 


DNApool 




R:AGCAGAGGCACACGAGATG 




partial IVSS; E6; partly IVS6 


(nt87924-88210) 


sequencing 


P7 


F:ACGTGTTTTCACTCTCAACCA 


60.0 


329 bp, including 


NC_010450.3 


DNApool 




TJTTP ATA AP A AP AnAnnnnnA 




pariiy i v ou, tj. / , paniai i v o / 




sequencing 


P8 


F:GCTGTTCTCCTCCCATCTCT 


60.0 


282 bp, including 


NC_010450.3 


DNApool 




R:ATGAAGCAGTTCCACGGTT 




partial IVS7; E8; 


(nt91524-91806) 


sequencing 








partial IVSS 






P9 


F:TCCCCTTTCTTTTTCAGTGG 


63.0 


234 bp, including 


NC_010450.3 


DNApool 




R:TCCCAGCATTAGATACACCCT 




partial IVSS; E9; partial IVS9 


(nt98070-98304) 


sequencing 


PIO 


F:TCTTGGTTGACTGCCTGCC 


59.0 


359 bp, including 


NC_010450.3 


DNApool 




R:CAGGGTTTGGGAAAGGAGC 




partial IVS9; ElO; partly IVSIO 


(ntl02365-10272) 


sequencing 


PIO 


F:AGCACGGAAGGAAAGACAGC 


58.0 


116 bp, including 


NC_010450.3 


DNA pool 




RiCAGCATTAGATACACCCTGGAATT 




partial IVS9; ElO; 


(ntl02489-10260) 


sequencing; 












PCR-RFLP 


Pll 


F:CATCAGTTGGTGAGGTCCTATT 


61.0 


322 bp, including 


NC_010450.3 


DNApool 




R:GAGCAAGAAGAGTCACAAGCA 




partial IVSIO; Ell; partial IVSll 


(ntl02979-11301) 


sequencing 


P12 


F:GACCCCAAAGGAGACTACC 


56.0 


813bp including 


NC_010450.3 


DNApool 




R:TTACAACAAGAGGGTGAGA 




partial IVSll; E12; 


(ntll7632-18445) 


sequencing 



PCR, polymerase chain reaction; LEF-1, lymphoid enhancer binding factor 1; Acc. No., means accession number; E, means exon; TVS, means intron; 



UTR, means untranscriptional region. 



sequenced by Invitrogen Corporation and then analyzed by 
the Chromas software. 

Polymerase chain reaction restriction fragment length 
polymorphism (PCR-RFLP) 

The PCR-RFLP technique was used to detect the 
polymorphism of 99514A>G and 119846C>T loci in 
porcine LEF-1 gene. Ahquots of 20 |iL PCR products were 
digested with 10 U EcoR Hand Hinfl (Takara) for 4 h at 
37°C following supplier instructions, respectively. The 
digested products were detected by electrophoresis on 2.0% 
agarose gel stained with EB. 

Statistical analysis 

Statistical analysis was performed on the basis of 
records of the number of teats in GZ, HJ, BM, and LW 
breeds. Genotypic and allelic frequencies of polymorphism 
loci on the LEF-1 gene were analyzed by the T-test for the 
four pig breeds using the SPSS software (version 19.0). 



Gene homozygosity (Hq), heterozygosity (Hg), effective 
allele numbers (Ng), and polymorphism information content 
(PIC) were evaluated by means of Nei's methods (Nei et al., 
1979). The association between genotypes and teat number 
trait was analyzed through the general linear model 
procedure of SPSS software (version 19.0) according to the 
following model: 

Yy = (i-HGenotypej-nEij 

Where Yy was the observation of animal, \i was the 
overall mean, genotypej was the fixed effect of the genotype, 
and Eij was the random effect. Effects associated with some 
of the factors, were not included in the above mentioned 
models, because the prehminary statistical analysis 
indicated that these factors were not significant sources of 
variation in the populations. Significant difference was 
stated when the genotype effect was a significant source of 
variation, and the p-value for the difference between the 
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least square means for each genotype was less than p = 0.05 
level. 

RESULTS 

Expression of lymphoid enhancer binding factor 1 in 
porcine different tissues 

In 3-day-old Guanzhong Black pigs, we detected the 
relative mRN A levels of porcine LEF-1 in teat, heart, liver, 
subcutaneous fat, spleen, kidney, lung and longissimus 
muscle tissues (Figure lA). The data showed a wide range 
of expression of LEF-1 in different porcine tissues. In detail, 
relatively high levels in teat, spleen, kidney and lung tissue 
(about 25.6, 7.3, 10.3, and 11.9 fold to heart tissue 
respectively) were observed As well as low levels in 
subcutaneous fat tissue and longissimus muscle tissue. 

In 180-day -old Guanzhong Black pigs, the relative 
mRN A levels of porcine LEF-1 in teat, heart, liver, spleen, 
lymph, thyroid and fat tissues are shown in Figure IB. High 
abundance was detected in teat tissue and low in heart tissue. 
In the important glands, such as thyroid, LEF-1 mRNA 
level was some 8 fold higher than heart tissue. In addition. 



spleen and liver also had a relatively high level (about 3 
fold to heart), while low in subcutaneous fat and lymph 
node. 

Furthermore, we detected the expression of LEF-1 in 
teat tissues of 3-day-old, 45-day-old and 180-day-old 
Guanzhong Black pigs as well (Figure IC), and the data 
displayed a decrease of the relative mRNA levels of porcine 
LEF-1 with age, in detail, the porcine LEF-1 was expressed 
at a level 3.5-fold greater in 3-day-old GZ pigs and 2.1-fold 
greater in 45-day-old GZ pigs than that in 180-day-old GZ 
pigs. 

Two novel genetic variations of porcine lymphoid 
enhancer binding factor 1 gene 

The porcine LEF-1 gene is located in SSC8 and 
contains 153883 bases including 12 exons and 11 introns 
(Figure 2). The LEF-1 mRNA consists of 2357 nucleotides 
and codes 398 amino acids. Two novel mutations 
(NC_010450.3: g. 99514A>G 119846C>T) were firstly 
detected by DNA pool sequencing in exon 4 and intron 9, 
respectively (Figure 3). No polymorphism was detected in 
the remaining sequences (P1-P3, P5-P9, Pll, and P12) 
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Figure 1. LEF-1 mRNA levels in Guanzhong Black pigs. (A) LEF-1 mRNA levels in 8 tissues of 3-day-old Guanzhong Black pigs, 
including teat, heart, liver, subcutaneous fat, spleen, kidney, lung and muscle tissues were detected. All the mRNA levels were 
normalized to LEF-1 mRNA of heart; (B) LEF-1 mRNA levels in 7 tissues of 180-day-old Guanzhong Black pigs, including teat, heart, 
liver, spleen, lymph, thyroid and fat tissues, were detected. All the mRNA levels were normalized to LEF-1 mRNA of heart; (C) 
Comparative analysis of LEF-1 mRNA expression in Guanzhong Black pigs teat tissues during 3 growth stages. All the mRNA levels 
were normalized to LEF-1 mRNA in 180-day-old. LEF-1, lymphoid enhancer binding factor 1. 
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Figure 2. The structure of porcine LEF-1 gene. Black boxes correspond to CDs. PI to P12 denotes the fragments amphfied by primers 
shown in Table 1. The SNP sites detected in present study were also included in this figure. LEF-1, lymphoid enhancer binding factor 1; 
CDs, coding sequence region; SNP, single nucleotide polymorphism. 



(Table 1) of the LEF-1 gene. The 99514A>G mutation Population genetic indexes 

resulted in a synonymous genetic code of proline, in detail. The population genetic indexes involving in Hq, Hg, Ng, 

CCA (Pro)>CCG (Pro) at position 166 aa of LEF-1. The and PIC are listed in Table 2. Generally, when PIC value 

119846C>T locus was observed in intron 9 and had no <0.25, it was defined as the low polymorphism; when 0.25 

effect on amino acids changing. <PIC value <0.5, it was median polymorphism; when PIC 

TCTCC4(G>GGiTC g 




Figure 3. Electrophoresis patterns of the DNA sequencing map and EcoR\l and Hinfl PCR-RFLP at porcine LEF-1 gene. (A) DNA 
sequencing map (upper) and electrophoresis pattern of the EcoR II PCR-RFLP (lower). The genotype "GG" demonstrated one undigested 
band (344 bp); Genotype "AA" showed two digested fragments (155 bp and 189 bp); Genotype "AG" showed three fragments (155, 189, 
and 344 bp). (B) DNA sequencing maps (upper) and electrophoresis pattern of the Hinf I PCR-RFLP (lower). The genotype "CC" 
demonstrated one undigested band (116 bp); genotype "TT" showed one digested fragment (89 bp) and; genotype "CT" showed two 
fragments (89 and 116 bp). Another digested fragment (27 bp) was invisible due to the lack of resolution within agarose gel. PCR-RFLP, 
polymerase chain reaction-restriction fragment length polymorphism; LEF-1, lymphoid enhancer binding factor 1. 
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Table 2. Genetic indexes of the 99514A>G and 119846C>T 
mutation within porcine LEF-1 gene 



Loci 


Breeds 


Ho 


He 


Ne 


PIC 


99514A>G 


GZ 


0.3787 


0.6212 


2.6406 


0.3537 




HJ 


0.8824 


0.1176 


1.1724 


0.1479 




LW 


0.6970 


0.3030 


1.3912 


0.2532 




BM 


0.3226 


0.6771 


3.0998 


0.3415 


119846C>T 


GZ 


0.4308 


0.5692 


2.3213 


0.3447 




HJ 


0.8387 


0.1613 


1.1923 


0.1860 




LW 


0.6875 


0.3125 


1.4546 


0.2289 




BM 


0.3000 


0.7000 


3.3333 


0.3747 



LEF-1, lymphoid enhancer binding factor 1; Ho, gene homozygosity; He, 
gene heterozygosity; Ne, effective allele numbers; PIC, polymorphism 
information content; GZ, Guanzhong Black pigs; HJ, Hanjiang Black 
pigs; LW, l^ge White pigs; BM, Bamei pigs. 



value >0.5, it was high polymorphism (Mateescu et al., 
2005). All breeds analyzed were in accordance with the 
Hardy- Weinberg equilibrium (p>0.05). In the 99514A>G 
mutation, low level PIC was observed in HJ and median 
level in GZ, LW, and BM; In the 119846C>T mutation, low 
level PIC was observed in HJ and LW and median level in 
GZ and BM; Additionally, the homozygosity occupied a 
higher rate in HJ and LW than that in GZ and BM. 

Genotypic and allelic frequencies 

The frequencies of widespread allele "A" in 99514A>G 
mutation ranged from 0.5536 to 0.9118, and "T" in 
119846C>T mutation ranged from 0.4833 to 08824 in four 
pig breeds (Table 3). Generally, more individuals with 
genotype "AA" and "TT" were identified in HJ and LW, yet 
few individuals with genotype "GG" emd "CC" were 
detected in GZ, HJ, LW, and BM. 

Association analysis of genetic variations with teat 
number trait 

The association analysis between different genotypes 
and teat number trait was assessed in GZ, HJ, BM, and LW. 



The 99514A>G mutation revealed an extremely 
significant statistical relevance between different genotypes 
and teat number trait in GZ (p<0.001) and LW (p = 0.002), 
and significant relevance in GZ (p = 0.017) (Table 4). On 
the whole, the individuals with genotype "GG" or "AG" 
showed more teat number than those with genotype "AA". 

The 119846C>T mutation suggested significant 
association of different genotypes with teat number trait in 
GZ (p = 0.042) and LW (p = 0.003) and no significance 
observed in HJ (p = 0.252) and BM (p = 0.591) (Table 4). 
In general, we fotmd individuals with genotype "TC" had a 
higher teat number than those with "TT" in the four breeds. 

Analysis of haplotype frequencies in Guanzhong Black 
and Large White pigs. 

To fiirther observe the correlation between the 
99514A>G and the 119846C>T mutation, the haplotype 
frequencies were analyzed in GZ and LW. The data showed 
that haplotypes "AGTC" and "AATT" accounted for a large 
proportion in GZ (34.85% and 57.58%) and LW (58.59% 
and 37.38%), while other haplotypes accounted for a 
proportion of less than 8% (Table 5). In addition, we found 
only one haplotype "GGCC" in LW and the individual had 
a higher porcine teat number than most of others with 
haplotypes "AGTC" and "AATT" (The data were not 
shown), the results implied that "G" and "C" could encode 
for a high teat number in pigs. 

DISCUSSION 

Here, the levels of LEF-1 mRNA were detected by 
qPCR in multiple tissues of 3 -day-old and 180-day-old GZ 
and the results both showed relatively high abundance in 
teat tissue. In addition, by comparing the expression of 
porcine LEF-1 among 3-day-old, 45-day-old, and 180-day- 
old Guanzhong Black pigs in teat tissues, we found the 
relative levels of porcine LEF-1 mRNA decreased with age. 



Table 3. Genotypic and allelic frequencies of 99514A>G and 119846C>T mutation within porcine LEF-1 gene 



Ix)ci 


Breeds 


Total - 




Genotypic frequencies 




Allelic frequencies 




AA(Paa) 


AG(Pag) 


GG(Pgg) 


A(Pa) 


G(Pg) 


9951 4 A>G 


GZ 


198 


69(0.3485) 


123(0.6212) 


6(0.0303) 


0.6591 


0.3409 




HJ 


34 


29(0.8529) 


4(0.1176) 


1(0.0294) 


0.9118 


0.0882 




LW 


99 


66 (0.6667) 


30 (0.3030) 


3(0.0303) 


0.8182 


0.1818 




BM 


30 


7(0.2333) 


19(0.6333) 


4(0.1333) 


0.5536 


0.4464 




Total 


361 


171(0.4736) 


176(0.4875) 


14(0.0388) 


0.7175 


0.2825 








TT(P„) 


TC(Ptc) 


CC(Pcc) 


T(Pt) 


C(Pc) 


119846C>T 


GZ 


198 


78(0.3939) 


111(0.5606) 


9(0.0455) 


0.6742 


0.3258 




HJ 


34 


26(0.7647) 


8(0.2353) 


0(0) 


0.8824 


0.1176 




LW 


99 


69(0.6970) 


30(0.3030) 


0(0) 


0.8485 


0.1515 




BM 


30 


4(0.1333) 


21(0.7000) 


5(0.1667) 


0.4833 


0.5167 




Total 


361 


177(0.4847) 


170(0.4663) 


14(0.0491) 


0.7258 


0.2742 



LEF-1, lymphoid enhancer binding factor 1; GZ, Guanzhong Black pigs; HJ, Hanjiang Black pigs; LW, Large White pigs; BM, Bamei pigs. 
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Table 4. Association of the polymorphic locus with teat number trait 



SNP 


Genotypes 




Teal numbers 


(meaniSD) 




GZ 


HJ 


LW 


BM 


99514A>G 


AA 


12.36+1.16 


13.52+0.87 


13.35±1.16 


11.81+0.77 




AG 


13.28+1.59 


14.20+1.10 


14.09+0.80 


12.29+1.18 




GG 


15.33±0.82 


16.00 


14.67±1.15 


13.00+1.55 




p value 


<0.001 


0.017 (<0.05) 


0.002 (<0.05) 


0.136 


119846C>T 


TT 


12.33+1.38 


13.48+0.92 


13.48±1.23 


12.33+1.74 




TC 


12.38+1.78 


14.33+2.89 


14.22±0.80 


13.25+1.71 




CC 


13.77±1.64 


0 


0 


12.60+0.89 




p value 


0.042 (<0.05) 


0.252 


0.003 (<0.05) 


0.591 



SNP, single nucleotide polymorphism; SD, standard deviation; GZ, Guanzliong Black pigs; HJ, Hanjiang Black pigs; LW, J^ge White pigs; BM, Bamei 
pigs. 



Collectively, it strongly suggested a regulatory potential in 
development of porcine breast. 

The development of porcine breast was mainly 
composed of the formation and speciaUzation of milk Une, 
formation of mammary placode, glandular bud and 
primitive ductal branching (Bahnsky, 1950). The above 
phases of breast development were complicated and tightly 
regulated by a series of signaling pathway including Wnt, 
fibroblast growth factor (FGF) and T-box 3 (Tbx3). 

In Wnt signaling pathway, several molecular markers, 
such as WntlOb and the Wnt reporter (TOP-Gal), play key 
roles in formation and specialization of mammary line. 
When dickkopf Wnt signaling pathway inhibitor 1 (Dkkl), 
inhibiting expression of TOP-Gal and WntlOb, was 
expressed, all morphologic evidence of mammary 
development disappeared (Chu et al., 2004; Veltmaat et al., 
2004; Veltmaat et al., 2006). 

The FGFIO is firstly expressed in the formative period 
of mammary line and regulates FGFR2b to affect Wnt 10 
response to the development of the mammary line. 
Additionally, FGFIO"'" and FGFR2b"'" mice showed 
developmental defects of mammary bud. At the same time, 
FGF/FGFRl signals were also required early in the 
initiation of mammary development (MaiUeux et al., 2002; 
Eblaghie et al., 2004). 

The Tbx3 is relevant to the ulnar-mammary syndrome 
in humans and haploinsufficiency of Tbx3 causes severe 
mammary hypoplasia, like the complete loss of mammary 
glands (Bamshad et al., 1997). Moreover, Tbx3"'" mice 

Table 5. Analysis ol'liaplolypc Ircqucncies in GZ and LW 



Breeds Haplotypes 



Numbers 



JJaplotype 
frequencies (%) 



GZ 


AATT 


69 


34.85 




AGTC 


114 


57.58 




Others 


3 


7.57 


LW 


AATT 


58 


58.59 




AGTC 


37 


37.38 




Others 


4 


4.03 



GZ, Guanziiong Black pigs; LW, Large White pigs. 



displayed an abnormal phenotype, which failed to develop 

most mammary placodes (Davenport et al., 2003; Jerome- 
Majewska et al., 2005). Thus, Tbx3 is required early in the 
initiation of mammeiry development. Furthermore, Wnt 
signaling and FGF/FGFRl signaling pathways converge to 
induce and maintain mesenchymal Tbx3 (Eblaghie et al., 

2004) . 

The LEF-1 gene could be treated as a bridge 
maintaining these three signaling pathways. In Wnt 
signaling pathway, LEF-1 gene is a downstream 
transcriptional mediator, which regulates breast 
development by regulating activity of Wnt signaling. In 
FGF signaling pathway, FGFR2b could affect expression of 
LEF-1 gene via regulating WntlO (MaiUeux et al., 2002; 
Eblaghie et al., 2004). Additionally, FGF signaling was at 
least in part mediated by LEF/TCF-P-catenin signaling 
(Holnthoner et al., 2002). In Tbx3 signaling pathway, 
preplacodal Tbx3 accumulation was essential for WntlOb 
and LEF-1 in development of mammary line (Davenport et 
al., 2003; Eblaghie et al., 2004; Jerome-Majewska et al., 

2005) and prodigious Tbx3 increased LEF-1 and WntlOb 
expression (Cho et al., 2006). Therefore, LEF-1 gene plays 
an important role in development of porcine breast and it is 
important to search for potentially polymorphic loci. 

At present, the LEF-1 gene variations associated with 
porcine teat number have not been thoroughly investigated. 
In previous study, there was a 1351T>C (NM_001129967.1) 
transition and a 1666A>C (NM_00 11 29967.1) trans version 
found in the 3' end of LEFl mRNA. The 1351T>C locus 
had highly significant associations with inverted teat trait 
and total teat number trait, while the 1666A>C mutation 
was significantly associated with total number of inverted 
teats in commercial populations (Tetzlaff et al., 2009). Here, 
we further explored potential polymorphic loci of LEF-1 
gene and found two novel single nucleotide polymorphism 
(NC_010450.3:g. 99514A>G 119846C>T) located in exon 
4 (P4) and intron 9 (P9), respectively. Then, through 
statistical analysis, we found that the 99514A>G mutation 
showed extremely significant relevance with total teat 
number trait in GZ (p< 0.001) and LW (p = 0.002), and 
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significant relevance in GZ (p = 0.017), but no significance 
in BM, while the 119846C>T mutation displayed 
significant association of different genotypes with teat 
number trait in GZ (p = 0.042) and LW (p = 0.003) and no 
significance observed in HJ (p = 0.252) and BM (p = 0.591). 
The individuals with genotype "AGTC" or "GGCC" had a 
greater teat number than those with "AATT". It insinuated 
the loci "G" and "C" could be associated with higher teat 
numbers. In addition, the 99514A>G mutation itself 
demonstrated a synonymous genetic code change of proline, 
in detail, CCA (Pro) >CCG (Pro) at position 166 aa of LEF- 
1 (398 aa), which might affect the efficiency of the 
translation changing phenotype to some extent. However, 
the precise mechanism of regulation of the porcine nipple 
numbers is still unknown and further study is urgent. 

According to analysis of PIC, the 99514A>G mutation 
within porcine LEF-1 gene displayed low PIC in HJ and 
median level in GZ, LW, and BM; the 119846C>T showed 
low PIC in HJ and LW and median level in GZ and BM. It 
indicated the difference of teat number trait from interbreed 
and intrabreed individuals, which could be conducive to 
assessing important polymorphic loci relevant with teat 
number trait. 

GZ, HJ, BM, and LW are important pig breeds in China. 
It is necessary to improve their multiple traits using 
conventional breeding methods combined with MAS. In 
recent years, teat number trait has been approved as an 
important indicator for litter size. The difference of numbers 
and developmental situations of porcine nipples might be 
caused by a major gene and several minor genes. However, 
the major gene affecting teat number trait was still unknown 
and many groups were attracted to focus on this study. In 
our study, we assessed two novel polymorphic loci of 
LEF-1 gene correlationed with teat number trait in GZ, HJ, 
BM, and LW. The data strongly suggested that the 
99514A>G mutation showed extremely significant 
statistical relevance with porcine teat number trait in GZ 
and LW and significant relevance in HJ, while the 
119846C>T mutation showed significant association in GZ 
and LW. It suggested that the LEF-1 gene polymorphism 
could be used as genetic markers to improve breeding for 
porcine teat number trait. 
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